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We present new constraints on the relativistic neutrino effective number N e g and on the Cos- 
mic Microwave Background power spectrum lensing amplitude A\. from the recent Planck 2013 
data release. Including observations of the CMB large angular scale polarization from the WMAP 
satellite, we obtain the bounds N e g = 3.71 ± 0.40 and A L = 1.25 ± 0.13 at 68% c.L. The Planck 
dataset alone is therefore suggesting the presence of a dark radiation component at 91.1% c.l. and 
hinting for a higher power spectrum lensing amplitude at 94.3% c.L. We discuss the agreement of 
these results with the previous constraints obtained from the Atacama Cosmology Telescope (ACT) 
and the South Pole Telescope (SPT). Considering the constraints on the cosmological parameters, 
we found a very good agreement with the previous WMAP+SPT analysis but a tension with the 
WMAP+ACT results, with the only exception of the lensing amplitude. 
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I. INTRODUCTION 

The recent precise measurements of the Cosmic 
Microwave Background (CMB hereafter) temperature 
anisotropies released by the Planck collaboration [T] are 
providing the tightest constraints on cosmological param- 
eters to date [2]. 

In this paper, we use this new dataset to constrain two 
parameters that affect the CMB "damping tail" regime, 
at small angular scales, namely the neutrino effective 
number N e g and the lensing amplitude A^, that, from 
previous experiments, have been reported as not consis- 
tent with the standard expectations [3J. 

We remind here that N e s effectively measures the num- 
ber of relativistic degrees of freedom at recombination 
and is related to the energy density in relativistic "dark" 
particles p v by: 



p v = 



N, 



cir 



Pi- 



(1) 



where p 7 is the CMB photon energy density, with value 
today p lfi ~ 4.8 x 1CT 34 gcm~ 3 . 

In the standard scenario, assuming three relativistic 
neutrino families, the expected value is iV e ff = 3.046. Ob- 
servation of a different value could point to new physics, 
related to the neutrino sector, such as non standard neu- 
trino decoupling, sterile neutrinos, etc., or to even more 
exotic physics, such as axions, extra dimensions, early 
dark energy (see e.g. [5], [1] and references therein). 

On the other hand, the A^ parameter is a phenomeno- 
logical parameter introduced in [6], that simply rescales 
the lensing potential: 



cr -> Ai.cn 



(2) 



where Cf is the power spectrum of the lensing field. 
The expected value for this parameter in the standard 
framework is A\. = 1. A value different from one could 



indicate either the presence of a systematic, or the pres- 
ence of new physics (see e.g. [5] and [7]). 

The previous CMB measurements obtained by the At- 
acama Cosmology Telescope (ACT, [5]) and the South 
Pole Telescope (SPT, [S]), when combined with the latest 
observations from the WMAP satellite (WMAP9, [TO]). 
have indeed provided different values for these two pa- 
rameters, that are in tension at the level of two stan- 
dard deviations. As showed in [3], the ACT dataset 
gives N cS = 2.85 ± 0.56 and A L = 1.64 ± 0.36 at 68% 
c.L, while the SPT dataset gives N eS = 3.72 ± 0.46 and 
A L = 0.85 ±0.13 at 68% c.L. 

Given this tension, it is certainly timely to investigate 
the constraints that can be obtained for N e g and A\. from 
the new data from the Planck satellite. 

The Planck collaboration has already presented results 
in [2] on N e g and A\. separately. Here we extend this 
analysis by varying 7V c g and A\. simultaneously, i.e. tak- 
ing into account the possible correlations between these 
two parameters as in [3] , and by properly comparing the 
results with the previous ACT and SPT measurements 
in the N e g-Ai, plane. 

Our paper is simply organized as follows: in the next 
section we describe the analysis method, in Section III 
we present our results also considering Baryon Acoustic 
Oscillation (BAO) surveys and H measurements, while 
in Section IV we derive our conclusions. 



II. DATA ANALYSIS METHOD 

Our main CMB dataset consists in the Planck public 
data release of March 2013 [TJ. We compare this dataset 
with the theoretical models using the CAMspec likelihood 
version 6.2 for high multipoles and the commander ver- 
sion 4.1 likelihood for low multipoles [TTj. We also con- 
sider the WMAP low-i likelihood for polarization [10] . 
This dataset is identical to the "PLANCK+WP" case 
presented in the Planck papers [H [TT] . 

For BAO surveys we include the following datasets: 
SDSS-DR7 [23 at redshift z = 0.35, SDSS-DR9 [13] at 



Parameter 


Planck+WP 


WMAP9±SPT 


WMAP9±ACT 


Q b h z 


0.02306 ±0.00051 


0.02264 ±0.00051 


0.02283 ± 0.00052 


fic/i 2 


0.1239 ±0.0054 


0.1232 ±0.0080 


0.110 ±0.010 


e 


1.04124 ±0.00077 


1.0415 ±0.0012 


1.0412 ±0.0025 


r 


0.095 ±0.015 


0.088 ±0.014 


0.090 ± 0.014 


n s 


0.996 ±0.018 


0.982 ±0.018 


0.969 ± 0.019 


log[10 10 A s ] 


3.111 ±0.034 


3.169 ±0.048 


3.174 ± 0.045 


iVcff 


3.71 ±0.40 


3.72 ± 0.46 


2.85 ±0.56 


A L 


1.25 ±0.13 


0.85 ±0.13 


1.64 ±0.36 


^A 


0.736 ± 0.022 


0.736 ± 0.023 


0.728 ± 0.025 


io[Gyr] 


13.08 ±0.38 


13.14 ±0.43 


13.90 ±0.55 


1 LjYl 


0.264 ±0.022 


0.264 ± 0.023 


0.272 ± 0.025 


H [km/s/Mpc] 


74.9 ±3.7 


74.6 ±3.7 


69.9 ±3.7 



TABLE I. Constraints at 68% confidence level on cosmological parameters from our analysis using Planck±WP, WMAP9±SPT 
and WMAP9±ACT. 
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FIG. 1. Comparison of the results for Planck±WP, WMAP9±SPT and WMAP9+ACT datasets in terms of the 1-D posterior 
distribution functions for the parameters -/V e ff (left) and A-l (right). 



z = 0.57 and WiggleZ [H] at z = 0.44, 0.60, and 0.73. 

Finally, we include the recent measurements for the 
Hubble constant H from the analysis of [13] and we refer 
to this dataset as HST. 

For the analysis method we use the publicly avail- 
able Monte Carlo Markov Chain package cosmomc [TB] 
which relies on a convergence diagnostic based on the 
Gelman and Rubin statistic. We use the latest version 
(March 2013) which includes the support for the Planck 
Likelihood Code vl.O (see |http : //c osmologist . info/ 
cosmomc/). The plots shown in this work are obtained 
via the python codes included in the cosmomc package. 

We run over the six-dimensional space of standard cos- 
mological parameters: the baryon and cold dark mat- 
ter densities fib and f2 c , the ratio of the sound horizon 



to the angular diameter distance at decoupling 9, the 
reionization optical depth r, the scalar spectral index 
ns, and the overall normalization of the spectrum As 
at k = 0.05 Mpc" 1 . We consider purely adiabatic initial 
conditions and we impose spatial flatness. In addiction to 
these parameters we let the number of neutrinos species 
(assumed massless) N e g and the lensing amplitude pa- 
rameter Al to vary, assuming the following flat priors: 
1.047 < A^off < 10 and 0.0 < A L < 4.0. 

In our runs, we also marginalise over the foreground 
parameters as in [51 [TT] . Since the correlations between 
the cosmological and foreground parameters is minimal, 
we do not report their values in this paper. The posteri- 
ors on foregrounds are in excellent agreement with those 
reported in [2]. 
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straints in the 2-D N e g — A^ plane. As we can see, the 
Planck+WP constraint on iV e fi is in impressive agree- 
ment with the previous WMAP9+SPT constraint and 
in tension with the WMAP9+ACT value. On the other 
hand, the constraint on A^ from Planck+WP is in bet- 
ter agreement with the WMAP9+ACT result but also 
consistent with the WMAP9+SPT constraint. 

It is also interesting to note that the Planck+WP 
dataset is in general in better agreement with the 
WMAP9+SPT dataset than the WMAP9+ACT on most 
of the parameters as the Hubble constant Hq, the matter 
density J7 m and the scalar spectral index ns- This could 
appear as in contradiction with the results presented in 
[5] where, on the contrary, a better agreement with the 
WMAP9+ACT dataset was reported. The reason of this 
different conclusion is due to the fact that while in [2] 
the comparison of Planck with ACT and SPT has been 
done by fixing A^ = 1 , here we let this parameter to vary 
freely. A similar conclusion has been reached by the SPT 
collaboration [T7] . 



FIG. 2. Comparison of the 2-D posterior distribution function 
from the Planck+WP, WMAP9+ACT and WMAP9+SPT 
datasets in the N e s — Al plane. 



III. RESULTS 
A. CMB data only 

We first discuss our results obtained using the 
Planck+WP dataset. Posteriors on the cosmological pa- 
rameters are shown in the first column of Table ^ 

As we can see, the Planck+WP dataset provides an 
indication for a larger value of both N e g and Al, with 
N cS = 3.71 ± 0.40 and A L = 1.25 ± 0.13 at 68% c.L. 

The constraints on N a g reported by the Planck collab- 
oration for the same dataset is N e g — 3.51 ± 0.39 at 68% 
c.L. However, this constraint is obtained with the con- 
dition Al = 1. The slightly larger value obtained in our 
analysis clearly shows that there is a small correlation be- 
tween these two parameters and that fixing A^ = 1 could 
slightly bias the constraints on N e g to smaller values. 

It is interesting to compare these results with the ones 
previously obtained in |3] for the ACT and SPT datasets, 
as shown in Table[l] As we can see, the Planck+WP result 
on N cS is perfectly consistent with the WMAP9+SPT 
constraint of N c g = 3.72 ± 0.46, while there is a tension 
with the WMAP9+SPT result on A L = 0.85+0.13. Vice- 
versa, the WMAP9+ACT constraint iV eff = 2.85 ± 0.56 
is clearly in tension with the Planck+WP result, while 
there is a better agreement with the bound on the lensing 
parameter A^ = 1.64 ± 0.36. 

In order to better demonstrate the consistency be- 
tween these datasets, we plot in Figure [l] the 1-D pos- 
terior distribution function for A^ and N e g coming from 
these three analyses, while in Figure [2] we report the con- 



B. Adding BAO and HST 

In this section we include the HST and BAO 
datasets. We report the results obtained from 

the Planck+WP+HST, Planck+WP+BAO and 
Planck+WP+HST+BAO analyses in the three columns 
of Table |Hj respectively. For the sake of brevity, we refer 
to the Planck+WP dataset as "CMB" . 

We show the 1-D posterior probability distributions 
for N e g and A^ in Figure [3] and the 2-D confidence re- 
gions for iV e ff and A^ in Figure [4] for all the dataset 
combinations discussed: CMB+HST, CMB+BAO and 
CMB+BAO+HST. 

Our results are in perfect agreement with those already 
presented in [2]. In fact, the introduction of the BAO 
dataset tends to drag down both the value of N e g and 
Ajj, i.e. to a better consistency with the standard ex- 
pectation; on the other hand, the inclusion of the HST 
dataset preserves the Planck+WP mean values of param- 
eters, reducing the error bars and therefore increasing the 
hints for new physics. In particular, for the CMB+HST 
case, both N e g and A^ assume larger values at more than 
95% confidence level. 

When both BAO and HST are combined, the final ef- 
fect is to lower the value of A^ = 1.17 + 0.10, which how- 
ever is still not in full agreement with the expected value 
of unity, and to lower the number of neutrinos species 
as well, giving N c g = 3.56 ± 0.27, which also remains at 
almost la away from the standard value. 



IV. CONCLUSIONS 

In this brief paper we have reported new joint con- 
straints on the neutrino effective number N c g and the 
CMB lensing amplitude A\, from the new Planck dataset. 



Parameter 


CMB+HST 


CMB+BAO 


CMB±BAO±HST 


n b h z 


0.022953 ±0.00035 


0.02246 ±0.00031 


0.02262 ± 0.00028 


Q c h 2 


0.1234 ±0.0050 


0.1232 ±0.0053 


0.1260 ±0.0049 


e 


1.04123 ±0.00077 


1.04112 ± 0.00078 


1.04085 ± 0.00075 


T 


0.094 ±0.014 


0.087 ±0.013 


0.089 ±0.013 


n 3 


0.992 ±0.011 


0.974 ±0.011 


0.9815 ±0.0088 


log[W w A 3 ] 


3.108 ±0.030 


3.093 ±0.030 


3.103 ±0.029 


N eS 


3.63 ±0.27 


3.35 ±0.31 


3.56 ±0.27 


A h 


1.24 ±0.12 


1.16±0.10 


1.17±0.10 


«A 


0.733 ±0.014 


0.706 ±0.011 


0.7119 ±0.0094 


to[Gyr] 


13.15 ±0.23 


13.47 ±0.28 


13.27 ±0.23 


^ i?7i 


0.267 ±0.014 


0.294 ±0.011 


0.2881 ± 0.0094 


fl"o[km/s/Mpc] 


74.0 ±2.0 


70.4 ±1.9 


71.8 ±1.6 



TABLE II. Constraints at 68% confidence level on cosmological parameters from our analysis using CMB+HST, CMB+BAO 
and CMB+BAO+HST. 
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FIG. 3. Comparison of the 1-D posterior distribution functions from the CMB-only (Planck+WP), CMB+HST, CMB+BAO 
and CMB+BAO+HST datasets for N cB (left) and A L (right). 



We have shown that the Planck+WP dataset is hint- 
ing for both a presence of dark radiation (at the level of 
91.1%) and for an anomalous amplitude for the lensing 
parameter (at the level of 94.3%). The Planck+WP con- 
straints on JVeff and other parameters, such as the Hub- 
ble constant and the matter density, are in very good 
agreement with those obtained from the WMAP9+SPT 
dataset. It is clearly worth to note that two very dif- 
ferent datasets provide an indication for a larger value of 
the effective neutrino number. In general, we found a ten- 
sion on the derived parameters between the Planck+WP 
and the WMAP9+ACT datasets. This clearly indicates 
that the inclusion of the ACT dataset in a combined 
Planck+ACT has to be carefully considered. 



is more consistent with the results obtained from the 
WMAP9+ACT dataset, which also provide a ~ 2a indi- 
cation for a larger value. However, since the same signal 
is not found in the trispectrum analysis |18) . the nature 
of this anomalous lensing amplitude needs further inves- 
tigation. Moreover, our analysis clearly demonstrates a 
correlation between A\ and the main cosmological pa- 
rameters. 

The hints for new physics from the Planck+WP 
dataset are confirmed when the HST measurements are 
included and are weakened when the BAO dataset is con- 
sidered. The CMB+HST+BAO analysis also suggests 
the presence of anomalous values but at smaller statisti- 
cal significance. 



The anomalous lensing amplitude from Planck+WP 



It will be probably duty of the next Planck data re- 
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lease, with the full mission and polarization data, to pro- 
vide more precise, CMB only, constraints on the neutrino 
number and the lensing amplitude and to confirm or fal- 
sify these current hints for new physics from Planck. 



Acknowledgements 



eff 



FIG. 4. Comparison of the 2-D posterior distribution func- 
tions from the CMB-only, CMB+HST, CMB+BAO and 
CMB+BAO+HST datasets in the Nag—Ah parameters plane. 
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